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Contract  number:  W81XWH-1 2-1 -0286 


Title:  Improving  the  Diagnostic  Specificity  of  CT  for  Early  Detection  of  Lung  Cancer:  4D  CT- 
Based  Pulmonary  Nodule  Elastometry 

Principal  Investigator:  Peter  G  Maxim,  PhD 

Introduction: 

In  this  project  we  are  addressing  a  shortcoming  of  existing  lung  cancer  screening  methods  by 
developing  a  CT  based  method  of  characterizing  a  mechanical  property  of  pulmonary  lesions, 
specifically  tissue  elasticity  (stiffness)  that  should  have  a  higher  specificity  than  purely  anatomic 
low-dose  CT.  It  is  the  aim  of  the  proposed  study  to  decrease  the  false  positive  rate  of  CT 
screening  by  analyzing  the  mechanical  properties  of  suspiciously  appearing  tissue  during  CT 
screening.  We  hypothesize  that  malignant  pulmonary  nodules  are  less  elastic  (stiffer)  than 
benign  nodules  and  that  this  difference  in  elasticity  can  be  used  to  differentiate  cancerous  from 
benign  nodules,  which  would  help  to  decrease  the  false  positive  rates  of  CT  screening.  A 
measure  of  elasticity  can  be  derived  from  high-resolution  4-dimensional  computed  tomography 
(4D  CT)  using  deformable  image  registration  algorithms.  Unlike  conventional  3D  CT  imaging 
that  results  in  a  static  image  of  the  scanned  anatomy,  4D  CT  incorporates  also  the  temporal 
changes  of  the  anatomy  caused  by  respiratory  motion,  yielding  a  CT  ‘movie’  that  allows  the 
evaluation  of  tumor  motion  and  the  calculation  of  the  elasticity. 

Body: 

Specific  Aim  1.  Deveiopment  of  deformabie  image  aigorithms  for  processing  the  4D  CT 
images  to  determine  the  eiasticity  of  maiignant  and  benign  puimonary  noduies.  (Dr. 
Maxim,  Tasks  1,  months  1-8) 

Task  1.  Deveiopment  of  the  software  for  deformabie  image  registration,  anaiysis  of 
the  DVF  and  the  catcuiation  of  the  eiasticity  parameter  (Matiab). 

The  software  will  be  developed  using  the  mathematical  package  Matiab  (The 
Mathworks  Inc.,  Natick,  MA).  Two  deformable  image  registration  algorithms  will  be 
used  (DIR''°'  and  a  method  based  on  optical  flow,  DIR°'").  The  resulting  displacement 
vector  fields  will  be  analyzed  and  an  elasticity  parameter  for  the  pulmonary  nodules  will 
be  calculated  (Dr.  Maxim,  months  1  -  8). 

Status  (Task  1): 

A  manuscript  describing  our  algorithm  and  its  validation  has  been  submitted  to  ‘Radiotherapy 
and  Oncology’  (Green  Journal)  and  is  attached  to  this  report. 
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Specific  Aim  2:  Validate  our  method  in  rat  models  of  human  lung  cancer  and  benign 
inflammatory  lesions.  (Dr.  Maxim,  Tasks  2-4,  months  3  -  24) 

Task  2.  Preliminary  experiments:  Establish  optimal  protocol  for  the  benign 

pulmonary  model  (granulomatous  inflammation)  and  study  growth  kinetics. 

2a.  Purchase  animals:  Rowett  rats,  A549  and  SK-MES-1  cells  from  American  Tissue 
Culture  Collection  (ATCC),  carbon  nanotubes  (catalogue  number  900-1 501,  lot 
GS1801),  SES  research  (Houston,  TX)  and  necessary  culturing  media.  (Dr.  Maxim, 
months  1-3) 

2b.  Inoculate  15rats  (Rowett  nude  rats)  with  carbon  nanotubes  and  follow  with  serial 
MicroCT  measurements  to  study  growth  kinetics  to  establish  the  time  for  nodule 
development  to  reach  desired  size.  (Dr.  Maxim,  15  rats  total,  months  3-6) 

Task  3.  Grow  orthotopic  model  of  lung  cancer  and  benign  lesions  and  follow  with 
serial  MicroCT  imaging:  preliminary  experiments  to  establish  protocol  and 
optimize  software 

3a.  Inoculate  10  rats  with  orthotopic  human  lung  cancer  cells  (A549,  left  lung)  and  carbon 
nanotubes  (right  lung)  (Dr.  Maxim,  months  7-9) 

3b.  Acquire  CT  images  at  peak-inhale  and  peak-exhale  using  a  small  animal  ventilator  (Dr. 
Maxim,  month  9-10) 

3c.  Analyze  CT  images  and  derive  elasticity  parameter  and  optimize  software  if 
necessary.  (Dr.  Maxim,  month  10) 

Task  4.  Grow  orthotopic  model  of  lung  cancer  and  benign  lesions  and  follow  with 
serial  MicroCT  imaging,  analyze  data 

4a.  Inoculate  remaining  40  rats  (A549  cells,  left  lung  in  Rowett  nude  rats)  and  follow  with 
CT  imaging  at  peak-inhale  and  peak-exhale  (Dr.  Maxim,  months11-13) 

4b.  Perform  simplified  analysis:  Delineate  malignant  and  benign  pulmonary  nodules 
and  measure  volumes  at  peak-inhale  and  peak-exhale.  Derive  elasticity  parameter 
based  on  the  ratio  of  the  volumes.  (Dr.  Maxim,  months  14-15) 

4c.  Analyze  acquired  CT  images  and  derive  elasticity  parameter  by  analyzing  the 

displacement  vector  fields  and  perform  statistical  analysis.  (Dr.  Maxim,  months  16-18) 

4d.  Repeat  experiments  and  analysis  with  second  cancer  cell  line  (SK-MES- 1 ),  50 
Rowett  rats,  (Dr.  Maxim,  months  18-23) 

4e.  Publish  animal  study  results  (Dr.  Maxim,  month  24) 
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Status  (Tasks  2,  3,  4):  Due  to  ongoing  repairs  and  upgrades  of  the  GE-MicroCT  scanner,  our 
proposed  experiments  are  being  delayed  again.  The  capability  of  acquiring  4DCT  images  is 
hampered  by  a  defect  in  the  GE-MicroCT  scanner.  Given  this  delay,  we  asked  the  DoD  for  a 
one-year  no-cost  extension,  which  was  recently  approved. 

We  were  able  to  to  successfully  generate  our  benign  model  using  the  talc  (instead  of  the 
proposed  carbon  nanotubes)  as  shown  in  Figure  1 . 


As  soon  as  the  capability  of  acquiring  4DCT  images  is  restored  on  the  GE-MicroCT  scanner,  we 
are  ready  to  acquire  and  analyze  the  mechanical  properties  of  the  tumors  and  the  benign  tissue 
in  accordance  with  the  proposed  method. 

Specific  Aim  3:  Validate  our  method  in  a  retrospective  review  of  over  200  4D  CT  scans 
from  patients  previously  treated  in  our  department.  (Dr.  Loo,  Task  5  months  1  -  20) 

Task  5.  Analyze  approximately  200  4D  CT  images  from  previously  treated  patients 

and  patients  recruited  within  the  funding  period. 

5a.  De-archive  all  previously  acquired  thoracic  4D  CT  scans  and  identify  suitable 
patients  for  the  study.  Our  institutional  data  (all  4D  CT  scans)  are  currently  stored  on 
DVD’s.  Data  will  be  de-archived  and  suitable  lung  cancer  patients  (patients  with 
benign  and  malignant  pulmonary  nodules)  will  be  identified.  (Dr.  Loo,  months  1-3) 
5b.  Identify  benign  and  malignant  pulmonary  nodules  to  be  included  in  the  analysis 
and  delineate  nodules  at  each  respiratory  phase.  (Dr.  Loo,  month  4) 

5c.  Perform  simplified  analysis  by  calculating  the  ratio  of  the  volumes  with  respect  to 
peak-inhale.  (Dr.  Loo,  months  5-8) 

5d.  Analyze  all  4D  CT  images  and  derive  elasticity  parameter  by  analyzing  the 
displacement  vector  fields  and  perform  statistical  analysis  (Dr.  Loo,  months  9-15) 

5e.  Analyze  data  from  new  patients  acquired  during  the  award  period  (Dr.  Loo, 
months  15-18). 

5f.  Publish  human  study  results  (Dr.  Loo,  months  19-20) 
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Status  (Task  5):  We  continue  to  de-archive  and  analyze  more  patients.  We  have  identified 
several  patients  with  benign  nodules  (in  addition  to  malignant  nodules  that  were  treated  in  our 
department).  From  those  patients  with  benign  nodules,  30%  showed  strong  motion  artifacts  in 
4DCT  that  aggravated  the  analysis.  23%  of  thus  far  analyzed  patients  have  very  small  benign 
nodules.  With  the  acquired  CT  resolution,  the  nodules  are  comprised  of  a  few  voxels. 

Deformable  image  registration  for  objects  with  such  few  voxels  is  inaccurate  and  ‘noisy’,  thus  the 
data  will  have  limited  value.  Dr.  Loo  has  started  with  the  delineation  of  the  benign  and  malignant 
nodules.  Data  will  be  processed  and  analyzed  shortly. 

Key  Research  Accomplishments: 

Our  first  aim  was  to  develop  and  validate  an  automated  software  package  for  determining  PN 
elasticity  against  a  manual  contouring  method,  and  preliminarily  assess  its  ability  to  distinguish 
malignant  tissue  by  comparing  the  elasticities  of  malignant  PN  with  those  of  the  lung.  This  work 
is  now  completed  and  a  manuscript  detailing  the  methodology  and  the  results  was  submitted  to 
‘Radiotherapy  and  Onclogy’  (and  included  to  this  report). 

Reportable  Outcomes: 

The  following  abstracts  have  been  selected  for  POSTER  presentation: 

1.  Mohammadreza  Negahdar,  Billy  W  Loo,  Maximilian  Diehn,  Lu  Tian,  Dominik 

Fleischmann,  and  Peter  G  Maxim,  “Automated  Tool  for  Determining  Pulmonary  Nodule 
Elasticity  to  Distinguish  Malignant  Nodules,”  ASTRO  2014 

A  manuscript  summarizing  our  initial  validation  was  submitted  to  Radiotherpay  and  Oncology. 
The  abstract  is  included  in  the  ‘Supporting  Documentation’  section. 

Conclusion: 

We  have  successfully  accomplished  specific  aim  1  of  the  proposed  study.  We  now  have  a 
functional  software  to  process  and  analyze  4DCT  images  to  distinguish  malignant  and  benign 
PN.  Despite  setbacks  in  time  because  of  upgrades  of  the  small  animal  equipment  ad  some 
difficulties  with  the  human  data,  we  believe  we  are  on  track  to  carry  out  the  proposed  research 
project. 
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Supporting  Data: 


Abstract  submitted  to  the  Annual  Conference  of  ASTRO  (2014): 


Automated  Tool  for  Determining  Pulmonary  Nodule  Elasticity  to  Distinguish  Malignant 

Nodules 

Purpose:  To  develop  and  validate  an  automated  method  of  determining  pulmonary  nodule  (PN) 
elasticity  against  a  manual  contouring  method,  and  preliminarily  assess  its  ability  to  distinguish 
malignant  tissue  by  comparing  the  elasticities  of  malignant  PNs  treated  with  stereotactic  ablative 
radiotherapy  (SABR)  with  those  of  the  lung. 

Methods:  We  analyzed  breath-hold  images  of  30  patients  with  malignant  PNs  who  underwent 
SABR  in  our  department.  A  parametric  nonrigid  transformation  model  based  on  multi-level  B- 
spline  guided  by  Sum  of  Squared  Differences  similarity  metric  was  applied  on  breath-hold 
images  to  determine  the  deformation  map.  The  Jacobian  of  the  calculated  deformation  map, 
which  is  directly  related  to  the  volume  changes  between  the  two  respiratory  phases,  was 
calculated.  Next,  elasticity  parameter  will  be  derived  by  calculating  the  ratio  of  the  Jacobian  of 
the  PN  to  the  Jacobian  of  a  1cm  region  of  lung  tissue  surrounding  the  tumor  (E-ROI)  as  well  as 
the  Jacobian  of  the  whole  lung  (E-Lung). 

Results:  For  the  first  group  of  15  patients  we  evaluated  the  volumetric  changes  of  PNs  and  the 
lung  from  the  maximum  exhale  phase  to  the  maximum  inhale  phase,  whereas  the  reverse  was 
done  for  the  second  group  of  15  patients.  For  the  first  group,  mean  and  standard  deviation  for  E- 
ROI  and  E-Lung  were  0.91  ±0.09  and  0.86±0.18,  respectively,  which  was  verified  by  the  manual 
method.  For  the  second  group,  E-ROI  and  E-Lung  were  1.34±0.27  and  1.57±0.51,  respectively. 
These  results  demonstrate  that  the  elasticity  of  the  PNs  was  less  than  that  of  the  surrounding 
lung  (p<0.0037). 

Conclusion:  We  developed  an  automated  tool  to  determine  the  elasticity  of  PNs  based  on 
deformable  image  registration  of  breath-hold  images.  The  tool  was  validated  against  manual 
contouring.  Preliminarily,  PN  elastometry  distinguishes  proven  malignant  PNs  from  normal  tissue 
of  lung,  suggesting  its  potential  utility  as  a  non-invasive  diagnostic  tool  to  differentiate  malignant 
from  benign  PN. 
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Cancer. 


Background  and  Purpose:  To  develop  a  noninvasive  method  for  determining  malignant  pulmonary  nodule 
25  (MPN)  elasticity,  validate  it  against  manual  contouring,  and  assess  its  ability  to  distinguish  MPN  from  normal  lung. 

Methods  and  Materials:  We  analyzed  breath-hold  images  at  extreme  tidal  volumes  of  23  patients  with  30  MPN 
treated  with  stereotactic  ablative  radiotherapy.  Deformable  image  registration  (DIR)  was  applied  to  the  breath-hold 
images  to  determine  the  volumes  of  the  MPNs  and  a  ring  of  surrounding  lung  tissue  (ring)  in  each  state.  MPNs  were 
also  manually  delineated  on  deep  inhale  and  exhale  images  by  two  observers.  Volumes  were  compared  between 
30  observers  and  DIR  by  dice  similarity.  Elasticity  was  defined  as  the  ratio  of  the  volume  change  of  MPN  compared  to 
ring. 

Results:  For  all  30  tumors  the  dice  coefficient  was  0.79±0.07  and  0.79±0.06  between  DIR  with  observers  1  and  2, 
respectively,  close  to  the  inter-observer  dice  value,  0.81±0.1.  The  elasticity  of  MPNs  was  1.24±0.26,  demonstrating 
that  volume  change  of  the  MPN  was  less  than  that  of  the  surrounding  lung. 

35  Conclusion:  We  developed  a  semi-automated  CT  elastometry  method  based  on  DIR  that  distinguishes  biopsy- 
proven  MPN  from  normal  lung.  Our  data  suggests  that  CT  elastometry  may  be  useful  in  distinguishing  malignant 


from  benign  nodules. 
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Introduction 


Elastometry  provides  tissue  characterization  that  may  help  distinguish  malignant  from  benign  tissues,  in  contrast  to 
static  morphological  imaging,  which  cannot  provide  this  degree  of  information  ^  Several  modalities  of  elasticity 
imaging,  mostly  based  on  ultrasound  have  been  proposed  and  applied  to  a  number  of  clinical  applications  including 
breast  carcinoma  metastatic  melanoma  head  and  neck  carcinoma  and  colorectal  carcinoma  However, 

45  since  air  is  opaque  in  ultrasound  images,  this  imaging  modality  has  been  considered  insufficient  for  adequately 
visualizing  lung  tissue.  Dynamic  CT  (including  cine-CT  and  four  dimensional  (4D)  CT)  images  the  lung  in  different 
states  of  expansion  that  could  be  used  to  derive  elasticity  information  in  contrast  to  static  three  dimensional  (3D) 
CT.  In  this  study,  we  propose  to  develop  pulmonary  nodule  elastometry  derived  from  CT  images  acquired  at  breath- 
hold  at  extreme  tidal  volumes. 

50  To  demonstrate  proof  of  principle  that  pulmonary  nodule  elastometry  can  be  determined  from  breath-hold  CT 
images  at  extreme  tidal  volumes,  we  analyzed  the  CT  images  of  23  patients  with  malignant  pulmonary  nodules 
(MPN)  previously  treated  with  stereotactic  ablative  radiotherapy  (SABR)  at  our  institution  and  compared  the 
volume  change  of  the  treated  MPN  with  that  of  the  surrounding  lung. 


55  Methods  and  Materials 

Patient  Selection 

From  all  patients  treated  on  an  institutional  protocol  of  stereotactic  ablative  radiotherapy  (SABR)  for  lung  tumors 
between  November  2011  and  October  2012,  23  patients  had  both  inhale  and  exhale  breath-hold  CT  scans  performed 
at  the  time  of  simulation  available  for  review  and  are  included  in  this  analysis. 

60 

CT  Imaging 

Deep  inspiration  and  natural  expiration  breath-hold  CT  scans  for  these  patients  were  acquired  on  a  Discovery  ST 
PET/CT  Scanner  (General  Electric  Medical  Systems,  Waukesha,  WI),  using  the  following  acquisition  settings:  120 


KVp,  110-195  mAs,  1.25  mm  slice  thickness,  0.97-1.36  mm  pixel  size,  500-700  mm  display  field  of  view.  Images 
65  were  reconstructed  using  either  the  built-in  Bone  Plus  or  Soft  convolution  kernel. 

Contouring 

The  original  treatment  plans  were  delineated  on  either  exhale  or  deep  inhale  and  the  clinically  used  contours  were 
taken  as  reference.  Two  observers  manually  delineated  the  MPN  on  the  opposite  respiratory  phase  using  our 
70  treatment  planning  software  (Eclipse  VI 1,  Varian  Medical  Systems,  Inc.,  Palo  Alto,  CA).  To  address  inter-observer 
variability,  the  MPNs  delineated  by  the  two  observers  were  compared  to  each  other  by  calculating  the  dice  similarity 
coefficient  voxel-wise 

Deformable  Image  Registration 

75  An  intensity-based  free-form  deformable  image  registration  (DIR)  workflow  was  developed  in  MIM  Maestro  (MIM 
Software  Inc.,  OH)  and  applied  on  extreme  tidal  volumes  to  determine  the  deformation  map  between  deep  inhale 
and  natural  exhale  images.  The  underlying  algorithm  parameters  were  inherent  in  the  program.  The  similarity  and 
smoothness  criteria  were  combined  into  one  energy  function,  which  was  minimized  in  the  registration  process 
Since  the  main  effect  of  respiration  on  MPN  is  the  transition  of  MPN  whereas  deformation  and  volume  change  of 
80  MPN  is  minimal  we  separate  effect  of  large  displacement  between  two  extreme  breath-hold  images  from  volume 
changes  of  MPN  itself  Therefore  workflow  starts  with  a  tumor  to  tumor  rigid  registration  with  manual  adjustment, 
followed  by  the  DIR  method  which  deforms  the  inhale  image  to  exhale.  Rigid  registration  and  manual  adjustment 
yield  the  advantage  of  robustness  to  the  MPN  location. 

85  DIR  Validation 

We  adopted  dice  similarity  coefficient  to  assess  the  quality  of  overlap  between  the  manually  delineated  contour  and 
deformed  contour  As  shown  in  Figure  1  (a),  for  patients  with  MPNs  originally  delineated  on  inhale,  DIR  was 


used  to  create  an  exhale  contour.  The  deformed  contour  was  then  compared  to  the  two  manually  delineated  exhale 
contours  using  the  dice  similarity  metric.  Similarly  and  as  shown  in  Figure  1  (b),  for  patients  with  MPNs  originally 
90  delineated  on  exhale,  the  two  observers  manually  delineated  contours  on  deep  inhale.  Next,  DIR  was  used  to  create 
exhale  contours  from  the  manually  delineated  contours  on  deep  inhale.  These  deformed  contours  were  then 
compared  to  the  clinically  delineated  contour  by  dice. 


Quantitative  Analysis  of  the  Deformation  Map 

95  Because  of  the  need  for  one-to-one  correspondence  between  material  points  during  continuous  deformation,  the 
calculated  Jacobian  determinant  is  required  to  be  non-zero  The  degree  of  regional  lung  expansion  is  measured 
using  the  Jacobian  determinant  of  the  deformation  map  which  is  directly  related  to  specific  volume  change 


Jn 


Vn+^Vn 

Vn 


(1) 


where  is  the  volume  of  voxel  element  n,  and  AV^  is  its  change  in  a  different  respiratory  cycle.  A  value  of  one 
100  implies  no  volumetric  changes,  while  a  Jacobian  determinant  greater  than  one,  or  smaller  than  one  implies  local 
tissue  expansion  or  local  tissue  contraction,  respectively.  It  therefore  follows  that  for  a  continuous  deformation  to  be 
physically  possible,  the  Jacobian  determinant  must  be  greater  than  zero.  We  took  advantage  of  this  fact  to  evaluate 
the  accuracy  of  the  calculated  deformation  map. 


To  validate  the  volumetric  changes  as  determined  by  the  DIR  method,  we  used  the  manually  delineated  MPNs.  The 
105  volumetric  changes  were  expressed  as  the  ratio  of  the  volumes  in  both  phases. 
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where  cp  is  the  manually  measured  of  volume  ratio  of  MPN,  Vinhaie  the  volume  of  the  MPN  at  deep  inhale  and 
V Exhale  is  the  volume  of  the  MPN  at  natural  exhale. 


110  Elasticity 


We  derived  an  elasticity  parameter  defined  as  the  ratio  of  the  volumetric  change  of  the  MPN  to  the  volumetric 
change  of  a  1  cm  ring  of  lung  tissue  surrounding  the  MPN,  Figure  2: 


\Jmpn-M 

\jRing~'^\ 


(3) 


where  is  the  calculated  volumetric  ratio  (exhale/inhale)  of  MPN  and  Jping  is  the  calculated  volumetric  ratio 
115  (exhale/inhale)  of  a  1  cm  ring  around  the  MPN.  The  normalization  to  the  1  cm  ring  around  the  MPN,  was  introduced 
to  remove  the  effect  of  MPN  location  and  amount  of  motion  in  the  lung,  as  both  the  MPN  and  the  1cm  ring 
surrounding  lung  tissue  undergo  the  same  force.  Based  on  the  size  of  tumors,  choosing  1  cm  ring  helps  get  fairly 
equal  contribution  of  lung  tissue  and  MPN  into  the  calculated  elasticity  parameter.  The  calculated  elasticity  bigger 
than  one  (d^  >  1)  shows  that  the  volume  changes  of  MPN  is  more  than  that  of  the  ring  and  the  calculated  elasticity 
120  smaller  than  one  (d^  <  1)  shows  that  the  volume  changes  of  MPN  is  less  than  that  of  the  ring. 


Statistical  Analysis 

An  equivalence  test  was  performed  to  assess  that  the  manually  measured  volume  ratio  was  equivalent  to  the  DIR 
calculated  volume  ratio  of  MPN.  A  paired  t  test  was  performed  to  compare  the  volume  change  of  the  normal  tissue 
125  and  the  tumor  tissue  with  both  tissue  samples  coming  from  the  same  patient.  All  analyses  were  performed  using 
SAS  version  9.4  (SAS  Institute  Inc.,  Cary,  NC,  USA). 


Results 

Patient  and  Tumor  Characteristics 

130  We  have  included  23  lung  cancer  patients  with  non-small  cell  lung  cancer  (NSCLC)  that  were  treated  at  our 
institution.  The  patient  and  tumor  characteristics  are  given  in  Table  1. 


DIR  Validation 


For  all  patients,  the  dice  similarity  coefficient  between  the  delineated  tumor  by  observer  1  and  that  of  observer  2  is 
135  0.81±0.1,  whereas  the  dice  value  is  0.79±0.07  and  0.79±0.06  between  DIR  and  observer  1  and  DIR  and  observer  2, 

respectively.  Figure  3  shows  an  example  of  both  delineated  and  deformed  contours  for  a  representative  patient  as 
well  as  the  distribution  of  the  calculated  dice  similarity  for  all  patients.  For  all  patients  we  calculated  the  Jacobian 
determinant  of  deformation  map  and  validated  the  calculated  volume  ratio  of  MPN  from  the  deep  inhale  phase  to  the 
natural  exhale  phase  against  the  manual  contouring.  The  calculated  Jacobian  determinant  of  deformation  map  within 
140  the  lung  is  always  greater  than  zero.  We  performed  an  equivalence  test  with  the  indifference  zone  defined  as  [-0.25, 
0.25].  The  mean  of  the  differences  of  the  paired  values  was  0.0293  with  a  95%  confidence  interval  of  [-0.0176, 
0.0761].  The  95%  confidence  interval  falls  completely  within  the  indifference  zone  thus  proving  equivalence 
(p<0.0001).  The  average  difference  between  the  manually  measured  volume  ratio  of  MPN  (cp)  and  calculated 
volume  ratio  value  via  DIR  was  0.02±0.15.  Figure  4  shows  the  calculated  volume  ratio  by  DIR  versus  manually 
145  measured  volume  ratio  of  MPN  for  all  patients  as  well  as  the  distribution  of  the  difference  between  the  calculated 
volume  ratio  and  manually  measured  volume  ratio  of  MPN  for  all  patients  which  mainly  concentrated  around  zero. 

MPN  Elasticity 

For  all  patients  of  our  study,  the  mean  and  standard  deviation  of  the  calculated  volume  ratio  of  the  MPN  and  that  of 
150  the  1cm  ring  surrounding  lung  tissue  was  0.91±0.18  and  0.75±0.13,  respectively.  The  mean  and  standard  deviation 
of  the  difference  between  the  calculated  volume  ratio  of  the  MPN  with  that  of  1  cm  ring  surrounding  lung  tissue  was 
0.16±0.15  (p<0.0001).  Figure  5  (a)  shows  the  calculated  volume  ratio  of  MPN  and  1cm  ring  surrounding  lung  tissue 
for  all  patients.  It  mainly  distributed  below  the  unity  line  demonstrating  that  the  volume  ratio  of  the  MPNs  are  less 
than  that  of  the  1  cm  ring  surrounding  lung  tissue. 

155  Finally,  based  on  our  definition  of  the  elasticity  parameter  in  Eq.  3,  the  mean  and  standard  deviation  of  the 
calculated  elasticity  (d^)  of  the  MPN  was  0.78±1.1.  Figure  5  (b)  shows  the  distribution  of  the  calculated  elasticity 
for  all  patients. 

Table  2  (supplementary  data)  shows  tumor  size  in  inhale  and  exhale  as  well  as  the  calculated  volume  ratio  of  MPN 


by  DIR  and  the  calculated  elasticity. 


160 


Discussion 

We  have  designed  and  developed  a  method  to  determine  the  elasticity  of  MPNs  based  on  deformable  image 
registration  of  breath-hold  CT  images.  Various  deformable  image  registration  methods  and  mechanical 

indices  have  been  surveyed  to  calculate  the  most  accurate  and  robust  deformation  map  as  well  as  differentiation 
165  metric,  respectively.  We  developed  a  workflow  that  processed  breath-hold  images  at  extreme  tidal  volumes  and 
calculated  the  volume  change  and  the  elasticity  of  MPNs. 

We  have  demonstrated  robust  performance  of  the  DIR  method  such  that  the  difference  between  DIR  and  each  of  the 
observers  is  comparable  to  the  inter-observer  difference  as  assessed  by  Dice  similarity.  Our  method  can  measure  the 
volume  ratio  of  MPNs,  which  can  be  either  expansion  or  contraction.  The  accuracy  of  this  measurement  is 
170  demonstrated  by  the  small  difference  between  the  manually  measured  volume  ratios  to  that  of  the  DIR  method 
(0.02±0.15). 

Based  on  the  presented  results  in  Figure  5  (a),  0  <  jRtng  <  Jmpn  <  1  which  demonstrates  that  MPNs  have  a  lesser 
volume  change  than  the  1  cm  ring  surrounding  lung  tissue  during  exhalation.  Thus,  one  expects  =  1^^^^  ^1  <  1 

\jRing-M 

which  has  been  shown  in  Figure  5  (b),  =  0.78  +  1.1).  This  demonstrated  that  MPNs  have  different 

175  characteristics  in  terms  of  elasticity  compared  to  the  surrounding  lung  tissue  during  respiration. 

Different  parts  of  the  lung  have  different  motion  amplitudes  and  volume  changes,  which  makes  the  calculated 
volume  ratio  of  MPNs  location-dependent.  However  the  formulated  elasticity  is  not  location-dependent  because  the 
normalization  of  the  volume  change  of  MPN  to  that  of  the  1  cm  ring  surrounding  lung  tissue  should  remove  the 
effect  of  MPN  location,  as  both  the  MPN  and  the  1cm  ring  surrounding  lung  tissue  undergo  the  same  force, 
180  irrespective  of  their  location.  Therefore,  it  allows  us  to  compare  the  elasticity  of  two  different  MPNs  in  two  different 
locations  of  lung. 

Our  study  does  have  certain  limitations.  First,  our  method  was  applied  to  a  small  number  of  patients,  and  focused 
only  on  proven  MPNs,  excluding  benign  nodules.  The  imaging  protocol  for  all  patients  was  not  completely  uniform 


which  resulted  in  different  image  resolutions  that  can  impact  the  analysis.  However  we  employed  image 
185  interpolation  so  that  DIR  was  applied  to  images  of  the  same  resolution. 

To  the  best  of  our  knowledge,  this  is  the  first  study  to  calculate  the  elasticity  of  MPN  from  dynamic  thoracic  images. 
This  provides  the  technical  framework  for  testing  the  hypothesis  that  MPN  and  benign  pulmonary  nodule  differ  in 
elasticity.  This  difference  may  come  from  increased  interstitial  fluid  pressure  in  tumors  which  can  change  the 
elasticity  of  MPN  prior  to  any  changes  in  morphological  features.  Studies  have  shown  that  most  solid  tumors  have 
190  increased  interstitial  fluid  pressure.  This  has  been  shown  for  breast  carcinoma  metastatic  melanoma  head  and 
neck  carcinoma  and  colorectal  carcinoma  with  values  as  high  as  60  mmHg.  The  tumor  interstitial  fluid  pressure 
is  uniform  throughout  the  center  of  the  tumor  and  drops  steeply  in  its  periphery.  The  mechanisms  that  determine 
the  increased  tumor  interstitial  fluid  pressure  are  not  fully  understood,  but  are  thought  to  involve  blood-vessel 
leakiness,  lymphatic  vessel  abnormalities,  interstitial  fibrosis  and  a  contraction  of  the  interstitial  space  mediated  by 
195  stromal  fibroblasts,  all  of  which  are  hallmarks  of  cancer 

Since  a  limitation  of  current  CT  based  screening  for  lung  cancer  is  distinguishing  malignant  from  benign  nodule 
our  method  could  impact  clinical  practice  by  increasing  the  specificity  of  CT-based  lung  cancer  screening.  In 
addition  CT  elastometry  may  predict  tumor  aggressiveness  that  could  potentially  be  used  in  multifocal  lung  cancer 
patients  to  ascertain  and  treat  the  most  aggressive  lesions  first.  Ultimately  it  may  even  be  helpful  to  distinguish 
200  aggressiveness  of  different  tumor  regions  to  guide  radiation  therapy  dose  painting  or  adaptive  boost  design. 

Future  directions  include  applying  our  analysis  to  both  proven  malignant  and  benign  nodules,  and  ultimately  may 
differentiate  the  aggressiveness  of  tumors  or  tumor  regions  to  aid  radiation  therapy  target  selection  or  dose  painting. 
Strain  tensor  map  along  with  elasticity  can  measure  pulmonary  nodule  homogeneity  in  order  to  compare  benign  and 
malignant  nodules.  In  addition  we  will  attempt  to  extend  our  method  to  4D  CT  images  routinely  acquired  for 


205 


radiotherapy  simulation. 
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Table  1.  Patient  and  tumor  characteristics. 


Characteristic  Value  (%) 

All  Patients 

23 

Gender 

Male 

13  (57) 

Female 

10(43) 

Age,  y 

Median 

66 

Range 

35-94 

Tumors 

30 

Primary/Secondary 

18/12 

Radius  of  Tumors,  cm  (treatment  planning) 

Mean 

1.06 

Range 

0.3-2.1 

Volume  of  Tumors,  cm^  (treatment  planning) 

Mean 

4.9 

Range 

0.11-38.8 

Volume  of  Tumors  in  opposite  phases,  crr\  (Obs  1) 

Mean 

7.6 

Range 

0.13-39.62 

Volume  of  Tumors  in  opposite  phases,  crr\  (Obs  2) 

Mean 

7.58 

Range 

0.13-39.74 

Table  2.  Raw  data  of  all  tumors  from  inhale  to  e^ale  for  current  study.  Second  and  third  columns  show  the 


average  of  manually  delineated  MPN  by  two  observers  in  inhale  ande^ale,  respectively.  Fourth  column  shows  the 
calculated  volume  ratio  of  MPN  by  DIR.  Fifth  column  shows  the  calculated  elasticity  of  MPN  by  Equation  (3). 


Tumor 

Tumor  Size 
(cm^)  (Inhale) 
Manual 

Tumor  Size 
(cm^)  (Exhale) 
Manual 

Calculated  volume 

ratio  of  the  MPN 

DIR 

Elasticity  (d^) 
Equation  (3) 

1 

17.48 

17.19 

0.96 

0.22 

2 

12.08 

12.07 

1.15 

1.35 

3 

0.139 

0.125 

0.68 

0.87 

4 

0.997 

1 

0.52 

1.42 

5 

1.74 

1.7 

0.98 

0.14 

6 

3.75 

3.85 

1.08 

1.55 

7 

8.61 

8.57 

1.00 

0.01 

8 

7.03 

7.54 

1.16 

0.75 

9 

22.62 

24.4 

1.03 

0.08 

10 

11.46 

11.23 

1.02 

0.06 

11 

7.18 

7.06 

0.90 

0.43 

12 

9.63 

9.77 

1.00 

0.04 

13 

3.69 

3.95 

0.97 

0.12 

14 

0.497 

0.462 

0.56 

0.93 

15 

7.33 

7.33 

0.93 

0.68 

16 

4.99 

4.45 

0.91 

0.32 

17 

3.42 

2.33 

0.94 

0.18 

18 

39.61 

38.36 

0.99 

0.63 

19 

19.53 

17.3 

0.89 

0.53 

20 

3.74 

4.1 

1.18 

5.94 

21 

2.45 

2.47 

1.26 

1.93 

22 

9.04 

8.91 

0.95 

0.2 

23 

0.81 

0.759 

0.69 

0.82 

24 

3.84 

2.31 

0.90 

0.17 

25 

1.72 

1.22 

0.69 

0.68 

26 

12.35 

10.58 

0.93 

0.32 

27 

0.39 

0.38 

0.73 

0.81 

28 

0.901 

0.864 

0.60 

1.12 

29 

1.32 

1.1 

0.72 

0.94 

30 


12.34 


12.02 


0.96 


0.14 


Inhale 


Exhale 


Clinically  delineated  Contour  (Reference) 


Manual  Delineation  (Obsl) 


Manual  Delineation  (0bs2) 


(a) 


Manual  Delineation  (0bs2) 


(b) 


figure  1.  (a)  Patients  with  MPNs  originally  delineated  on  inhale;  DIR  was  used  to  ereate  an  e>diale  eontour.  The 
deformed  eontour  was  then  eompared  to  the  two  manually  delineated  exhale  eontours  using  the  diee  similarity 
metrie  and  average  of  the  diee  values  has  been  reported,  (b)  Patients  with  MPNs  originally  delineated  on  e>diale; 
DIR  was  used  to  ereate  two  exhale  eontours  from  the  manually  delineated  eontours  on  inhale.  Then  the  deformed 
eontours  were  eonpared  to  the  elinieally  delineated  eontour  by  diee  and  average  of  the  diee  values  has  been 


reported. 


Deep  Inhale  CT  Natural  Exhale  CT  with  Deformed  Contours 


Figure  2.  An  axial  view  of  a  representative  patient;  delineated  MPN  (red)  and  lem  ring  of  the  surrounding 
lung  (yellow)  for  deep  inhale  (left  eolumn)  and  deformed  MPN  and  ring  eontours  superinposed  on  exhale 
image  (right  eolumn)  for  the  purpose  of  elastieity  ealeulation. 


Figure  3.  (a)  A  coronal  view  of  a  manually  delineated  tumor  in  exhale  of  a  representative  patient;  delineated 
MPN  in  exhale  (green)  and  delineated  MPN  in  inhale  deformed  into  the  exhale  (magenta),  Dice  similarity  for 
this  patient  is  0.791.  (b)  Distribution  of  the  calculated  Dice  similarity  metric  between  the  manually  delineated 
MPN  by  first  observer  with  that  of  second  observer,  (c)  Distribution  of  Dice  similarity  between  the  delineated 
MPN  in  inhale  deformed  into  the  e>diale  phase  by  DIR  with  the  manually  delineated  MPN  by  first  observer  in 
exhale,  (d)  Distribution  of  Dice  similarity  between  the  delineated  MPN  in  inhale  deformed  into  the  exhale  phase 
by  DIR  with  the  manually  delineated  MPN  by  second  observer  in  exhale.  Based  on  dice  similarity,  we  conclude 
that  DIR  conpares  favorably  to  manual  contouring  within  inter-observer  differences. 
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figure  4.  (Up)  Calculated  volume  ratio  of  MPN  versus  manually  measure  of  volume  ratio  of  MPN  (Equation 
(2))  for  all  30  patients  from  inhale  to  exhale.  Q'een  dashed  line  shows  unity  Jacobian  value  equivalent  to  no 


volume  ehanges.  Red  line  represents  the  identity  line.  (Bottom)  Distribution  of  the  differenee  between  the 
ealculated  volume  ratio  by  DIR  from  inhale  to  exhale  with  manually  measured  volume  ratio  of  MPN 


(0.02±0.15). 


(a) 


Elasticity 

(b) 


Figure  5.  (a)  Calculated  volume  ratio  of  1  cm  ring  of  surrounding  lung  (0.75±0.13)  versus  that  of  tumor 


(0.91±0.18)  for  all  patients.  Q*een  dashed  line  shows  unity  Jacobian  value  equivalent  to  no  volume  changes.  Red 


line  represents  the  identity  line,  (b)  Distribution  of  the  calculated  elasticity  by  Equation  (3)  from  inhale  to  exhale 


(0.78±1.1). 


